AD-A118  346  AEROSPACE  CORP  £L  SE6UND0  CA  CHEMISTRY  AND  PHYSICS  LAB  F/6  11/6 

AN  ANALYTICAL  MODEL  FOR  RF  INDUCTIVE  GLOW  DISCHARGE  COLUMNS  IN  — ETC(U) 
JUL  82  H  U  ECKERT  F04701-81-C-0082 

UNCLASSIFIED  TR-0062 (2472-06) -1  SD-TR-82-47  NL 


SD-TR-82-47 


AD  A118346 


REPORT  SD-TR^^r 


■  J> 

An  Analytical  Model  rf  Inductive 
Glow  Discharge  Colum^  in  Metal-Gas 
Mixtures  and  Application  to  the 
Rb-Xe  Systckn 


HANS  U.  ECKERT,  Consultant 
Chemistry  and  Physks  Labcmtory 
Laboratoiy  Operations 
The  Aerospace  Corporation 
El  Seguhdo.  Calif.  90245 


12  July  1982 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNUMITID 


r 


I 


DTTC 

.  .M - 

AU61  71880 

H 


Prepared  for 

SPACE  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
Los  Angeles  Air  Force  Station 
P.O.  Box  92960,  Worldway  Postal  Center 

Los  Angeles,  Calif.  90009  ^  ^  ^ 

82  08  16 


192 


This  report  was  submitted  by  The  Aerospace  Corporation,  El 
Segundo,  CA  90245,  under  Contract  No.  F04701-81-00082  with  the  Space 
Division,  Deputy  for  Technology,  P.O.  Box  92960,  Uorldway  Postal  Center,  Los 
Angeles,  CA  90009.  It  was  reviewed  and  approved  for  The  Aerospace  Corporation 
by  S.  Feuerstein,  Director,  Chemistry  and  Physics  Laboratory.  Lt  Efren  V. 
Fornoles,  SD/YLXT,  was  the  project  officer  for  Technology. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PAS)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
report's  findings  or  conclusions.  It  Is  published  only  for  the  exchange  and 
stloulation  of  Ideas. 


Efpen  V.  Fornoles,  1st  Lt,  USAF 
Project  Officer 


FOR  THE  CCMMANDER 


Jimmlfi  H.  Butler,  Colonel,  USAF 
DireJItor  of  Space  Systems  Technology 


REPORT  DOCUMENTATION  PAGE 


SD-TR-82-47 


a.  oovT  AceiaMON  mo 


4*  TtTUC  firffffftj 

AN  ANALYTICAL  MODEL  FOR  RF  INDUCTIVE  GLOW 
DISCHARGE  COLUMNS  IN  METAL-GAS  MIXTURES  AND 
APFLICATION  TO  THE  Rb-Xe  SYSTEM 


authomto 


f.  TVMi  or  nbmowt  •  rcwoo  eovinco 


«.  rurORMINO  ONO.  RBMONT  NUMOBR 

TR-0082 (2472-06)-! 


Hans  U.  Eckert 


F0470I-81-C-0082 


•.  reiiroRMiNO  oroanization  mams  amo  adorsm 

The  Aerospace  Corporation 
El  Segundo,  California  90245 


10.  RROORAM  KLIMIMT.  RROJKCT.  TASK 
ARKA  A  WORK  UNIT  NUMaCRS 


1 1.  eoNTROLUNO  orrici  namc  amo  aoorsss 

Space  Division 

Air  Force  Systems  Command 

Los  Angeles,  California  90009 


IS.  RSRORT  DATS 

12  July  1982 


la.  NuiaasR  or  RAoes 

39 


I  cmmlllmt  OMm^  IS.  SSCURITV  CLASS,  (at  SW«  irrmO 

Unclassified 


IS.  nSTRIOUTlON  STAT 


Approved  for  public  release;  distribution  unlimited. 


rr.  CNSTRISUTION  STATIM«NT  »•  tktumt  te  marnt  M,  It  «Mh«M  *mi  IUrhO 


to.  KBV 


MMMliy  Sr  MMk 


Alkali  Vapor,  Alkali  Vapor  Noble  Gas  Mixture,  Electroless  Discharge  Lamp, 
Glow  Discharge  Coluan,  Rubidium-Xenon  Spectral  Lamp 


The  electron  balance  equation  is  solved  for  a  glow  discharge  column  excited  by 
rf  induction  in  a  metal  vapor  -  noble  gas  mixture  where  the  metal  atoms  are 
ionised  in  a  single  step  and  ionisation  of  the  gas  is  negligible.  The  accumu¬ 
lation  of  metal  atoms  near  the  wall  is  accounted  for  by  postulating  for  their 
radial  distribution  n(r) ,  a  power  law  with  positive  exponent.  Under  this  con¬ 
dition,  the  profile  for  the  electron  density  n,(r)  can  be  expressed  by  zero 
order  Bessel  functions  containing  a  profile  parameter  q  that  is  related  to  thp 
exponent  in  the  atomic  distribution^ACr) ,  oils devaluated  as  a  function  of  rf 
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power  Input  and  electron  temperature  for  a  fixed  coljjmn  radius  R  «  0.4  cm 
by  solving  the  electron  energy  balance  equation  for  l!g»*f(r).  The  profiles 
n(r)  and  n^Cr)  are  Interrelated  by  the  condition  that  In  a  j^ady  state  their 
slopes  at  r  >  R  must  be  proportional  to  the  amblpolar  and  atomic  diffusion 
coefficients,  respectively.  This  and  the  common  boundary  condition  n^(R)  >  0 
provide  additional  relationships  that  permit  presentation  of  q  as  a  function 
of  the  two  control  variables:  power  Input  per  unit  column  length,  F/L,  and 
mean  concentration  of  metal  atoms  and  positive  Ions,  n^,  (determined  by  the 
vapor  pressure).  The  equations  are  applied  to  a  Rb-Xe  mixture  with  p(Xe)  ■ 

2  Torr.  The  Ionization  rate  coefficient  for  Rb  Is  derived  approximately  from 
measurements  of  Ionization  cross  sections  assuming  a  Maxwellian  distribution  f 
for  the  electron  energies.  Calculations  of  n  (r)  and  n(r),  carried  out  for 
two  values  of  P/L,  Indicate  a  tendency  for  thS  profiles  to  assume  a  rectan¬ 
gular  shape  %rlth  Increasing  F/L,  which  agrees  qualitatively  with  observations 
In  Cs  -  Ar  dc  columns.  The  values  of  T  associated  with  the  calculated 
cases  are  higher  than  observed  In  Rb-Xe^rf-exclted  lamps  and  suggest  that 
multistep  Ionization  may  not  be  negligible. 
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I.  INTRODUCTION 


Dc  and  ac  diachargea  in  mixturea  of  a  metal  vapor  with  a 
noble  gaa  have  found  wide  applications  in  space  illumination  and  their 
properties  have  been  studied  extensively  by  experiment  and  theory  [1- 12]. 
More  recently,  low  power  versions  of  such  discharges  have  been  used 
for  optical  pumping  in  various  devices  including  atomic  frequency 
standards  [1 3- 19],  In  this  case  the  metal  is  usually  rubidixim,  and  it 
has  been  foiud  preferable  to  excite  the  discharge  electrodelessly  by 
rf  induction  [15,  16].  Neither  the  gaseous  electronic  properties  of 
this  metal  nor  the  details  of  this  excitation  method  are  well  known, 
so  that  the  understanding  of  diis  kind  of  discharge  is  in  a  rudimentary 
state  and  its  development  has  been  carried  out  largely  on  empirical 
grounds.  Ofei  the  other  hand,  there  are  high  demands  on  stability, 
efficiency  and  useful  life  for  these  discharge  lamps,  notably  in  space-borne 
eqtiipment,  which  necessitate  a  thorough  understanding  of  the  ongoing  pro¬ 
cesses.  In  this  paper  a  first  step  toward  such  understanding  is  made 
through  establishment  of  a  simple  matheznatical  model  that  correlates 
some  internal  parameters  of  the  discharge  with  external  parameters  and 
yields  closed  form  expressions  for  the  radial  distributions  of  electrons  and 
metal  atoms.  The  metal  vapor  complicates  the  situation  compared  to  that 
in  a  simple  gas  discharge  by  its  tendency  to  segregate  from  the  buffer  gas. 
This  is  due  to  the  fact  that  the  metal  atoms,  having  a  considerably  lower 
ionization  potential  than  the  gas  atoms,  provide  most  of  the  ionization, 
although  their  average  concentration  is  generally  lower  by  orders  of 
magnitude.  The  metal  ions  travel  with  the  fast  electrons  in  ambipolar 
diffusion  to  the  wall  where  they  recombine  and  return  as  atoms  into  the 
discharge  volume.  Since  this  back  diffusion  occurs  at  a  slower  rate, 
the  metal  atom  concentration  builds  up  near  the  wall  and  may  exceed  here 
their  average  concentration  several  times.  It  is  of  importance  to  know 
the  metal  atom  concentration  directly  at  the  wall,  since  a  small  fraction 
is  retained  there  and  lost  from  the  discharge  process  through. some 
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physical  or  chemical  process.  As  long  as  the  retained  atoms  can  be 
replaced  by  vaporization  from  a  pool,  their  number  in  the  volume 
remains  constant  and  the  radiation  output  steady.  Once,  however, 
this  supply  is  exhausted,  the  output  declines  and  the  lamp  fails.  Being 
able  to  correlate  the  lamp  life  with  the  wall  atom  concentration  would 
help  in  a  quantitative  description  of  the  interaction  process.  Since  this 
concentration  is  difficult,  if  not  imxx>ssible  to  measure,  one  depends  on 
calculations,  for  which  the  model  will  be  useful. 

Knowledge  of  temperature  and  particle  distributions  is  also 
important  for  understanding  the  radiation  characteristics  of  the  lamp, 
notably,  the  degree  of  self  absorption  experienced  by  the  resonance 
radiation  before  leaving  the  plasma. 

The  operational  stability  of  the  lamp  also  depends  on  the  dis¬ 
tribution  of  particles  and  electric  fields.  It  is  well  known  that  a  low- 
pressure  rf  induction  discharge  at  small  power  inputs  usually  starts 
in  a  capacitive  mode;  ujMn  an  increase  of  power  it  abruptly  switches 
to  the  inductive  mode  and  vice  versa.  The  presence  of  two  species 
with  different  ionization  potentials  enhances  the  possibility  for  mode 
switching. 

So  it  is  eiq)ected  that  the  model  developed  in  the  following, 
which  to  the  author's  knowledge  represents  the  first  closed-form 
I  treatment  for  this  kind  of  discharge,  can  help  on  a  variety  of  problems. 
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n.  ANALYSIS 


a)  Radial  Distributions  of  Electrons  and  Metal  Atoms 


We  consider  an  axially  uniform  plasma  column  of  radius  R 
that  contains  free  electrons,  metal  atoms  and  singly  charged  ions  at 
concentrations  n^,  n  and  n.,  respectively,  which  are  functions  of  the 
radial  coordinate  r.  Noble  gas  atoms  which  are  present  at  concentra¬ 
tions  higher  by  several  orders  of  magnitude  do  not  participate  in  the 
ionization,  but  determine  the  diffusion  rates  of  atoms  and  charged 
particles.  As  mentioned,  there  is  a  counterflow  between  electrons 
and  ions  moving  to  the  wall  and  recombined  atoms  returning  into  the 
plasma.  Since  any  new  electron-ion  pair  generated  means  the  loss 
of  an  atom,  the  divergences  of  the  two  flows  are  oppositely  equal 

t^i^e  =  =  -Vr  {1),(2) 


where  j/^  is  the  ionization  frequency  of  an  average  electron  and  Fg  and  F 
are  the  densities  of  electron  and  atom  fluxes,  respectively.  Since  both 
fluxes  start  with  zero  at  their  respective  origins  r  =  0  and  r  s  R,  we 
also  have 


r*  =  -r 


(3) 


For  the  one  dimensional  column 


and. 


dne 

(4) 

e  “ 

d  r 

dn 

(5) 

i  w  -  JJ 

dr  • 

Combining  Eqs.  (4)  and  (5)  with  Eq.  (3)  yields 


(6) 


a  result  obtained  already  by  Waszink  and  Polman  [4].  D  and  Da> 
the  atomic  and  ambipolar  diffusion  coefficients,  which  are  essentially 
functions  of  the  gas  temperature  T  and  the  electron  temperature  T^, 
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are  considered  to  be  independent  of  radial  position.  With  Eq.  (4)  we 
can  write  the  electron  balance  equation  (1) 


d^n  ,  dn  y. 

+  4^”  =  0 

j  2  r  dr  Da  e 

dr 


is  proportional  to  the  local  value  of  n  and  can  be  expressed  by 


"i  =  ''gi“ 

where  is  the  ionization  rate  coefficient  of  ground  state  atoms.  It 
is  a  function  of  T<,.  The  distribution  of  n  is  assumed  to  follow  the 


power  law 


n  =  n.. 


r 


n) 


where  n_  is  the  value  of  n  at  r  =  R  and  q  >0.  Using  Eq.  (9)  and 

X\ 

setting  r/R  ■  p  we  can  write  Eq.  (7) 

^  °e  .  1  ^°e  .  „2  *'iRp  2q  =0 

j  p  2  p  d  P  Da 


where 


The  solution  of  Eq.  (10)  can  be  written 


.1/2 


T  r  R  /_1R  \ 
^e  ■  “eo  o  [q+1  \  Da  / 


where  n.  is  the  value  of  n  at  r  =  0  and  J  is  the  Bessel  ftinction  of 
®o  e  o 

zero  order.  Setting  at  the  wall,  P=  1,  n^  =  0  which  is  admissible,  if 
the  electron  mean  free  path  f  <<R,  we  obtain  Jq(X)  =  0  where 


X=  2.405  = 


R  /  *^^R  \  ^  ‘ 
1+1  \  Da  / 


t, 


Besides  at  the  axis,  where  it  is  automatically  satisfied,  condition  (6) 
will  be  satisfied  at  the  wall.  This  yields  by  differentiation  of  £qs.  (9) 
and  ( 1 1 )  fo  r  P  -  I 


2q  Hr  =  n^_  (q+1) 


where  Jj  is  the  first-order  Bessel  fvinction.  Replacing  nR  by  the  average 
atom  density  n,  which  through  Eq.  (9)  is  related  to  Ur  by 


Hr  =  n  (q  +1) 


Eq.  (13)  reduces  to 


^9“  =  “eo  “  ^  (X)  (IS 

Likewise  the  peak  electron  density  n^^  can  be  replaced  by  the  mean  value 


where 


“e  =  2neoIl(9) 


C]  (q)a/jo(XP9+ljpdP 


is  plotted  in  Fig.  1. 


Because  of  charge  neutrality  Qg  = 


we  have 


”e  =  ”t 


where  n^  is  the  total  average  concentration  of  metal  atoms  and  ions. 
It  is  given  by  _  ^ 


P  (T  .  ) 
*^v'  min 
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where  p  (T  .  )  is  the  metal  vapor  pressure  at  the  point  of  minimtim 
V  nun 

temperature  on  the  inside  wall  surface  and  k  is  the  Boltzmann  constant, 
n^  is  considered  to  be  one  of  the  input  variables.  The  other  independent 
variable  is  the  power  dissipated  in  the  discharge  column,  for  which  an 
expression  will  be  derived  in  the  next  section. 

b)  Dissipation  of  rf  Power 

The  input  of  power  P  into  the  discharge  col\imn  occurs  inductively 
by  an  rf  magnetic  field  of  angular  frequency  v  The  power  dissipated 
per  unit  of  coluzxui  length  L  follows  from  the  common  expression 


=  ZTTR^faE^pdp 


where  the  electrical  conductivity  <7  and  the  induced  electric  field  E  are 
both  ftinctions  of  p.  In  an  rf  field  the  real  part  of  <T  is  related  to  n^,  u 
and  the  electron-neutral  collision  frequency  by 


+  u^) 


where  e  and  m^  are  the  electronic  charge  and  mass,  respectively. 
The  variation  of  E  can  be  expressed  to  good  approximation  [20j  by 


E  =  Ej,pP 


where  E^^  is  the  value  of  E  at  the  boundary. 


(22) 


The  exponent  p  is  a  fiinction  of  the  ratio  R/  d,  where  d  is  the  skin  depth. 
For  R/6  <1,  pasl;for  R/^>1,  p  — >  R/d  .  d  is  defined  by 


(5  = 


1/2 


where  n  is  the  magnetic  permeability  of  the  vacuum. 

^  n 


(23) 


Although  E  has  a  strongly  nonuniform  distribution  across  the  column, 
it  has  been  found  by  several  experimenters  [21-23]  that  the  electron  temp¬ 
erature  in  low-pressure  induction  discharges  with  noble  gases  is 
little  affected  by  this  nonuniformity  and  can,  to  a  good  approximation, 
be  taken  as  constant.  There  may  be  two  different  reasons  for  this  result. 
In  the  first  place  it  is  apparently  a  consequence  of  the  high  thermal  con¬ 
ductivity  of  the  electron  component  which  averages  out  nonuniformities 
in  T^.  Secondly,  there  is  usually  a  capacitive  electric  field  from  the 
inductor  coil  present  in  the  discharge  which,  unlike  the  induced  field, 
does  not  vanish  at  the  axis.  Although  its  contribution  to  the  total  power 
input  is  small,  it  has  a  marked  effect  on  the  equalization  of  T^,  as  was 
noticed  by  Turban  [Z4]  for  a  discharge  in  hydrogen. 


We  thus  have  a  justification  for  combining  the  electron  density 
profile  Eq.  (11),  obtained  for  uniform  T^,  with  the  nonuniform  field 
distribution  Eq.  (22)  into  the  expression  for  the  power  input  Eq.  (20). 
Making  also  use  of  Eqs.  (11),  (12)  and  (21),  we  obtain 


1 


(24> 


In  many  practical  cases  Eq.  (24)  can  be  simplified.  Since  the  vapor 

pressure  of  the  metal  is  usually  in  the  milli-Torr  region  and  by  several 

orders  of  magnitude  lower  than  that  of  the  buffer  gas,  the  effective  degree 

of  ionization  is  small,  even  if  the  vapor  is  fully  ionized.  For  values  of 

n  in  the  order  of  10^2  _  lo^^cm"^,  y  10^  -  10^0  sec"^,  a»  —  10®-10^  sec 
e  c 


10 


and  R  1  cm,  it  follows  from  what  is  said  above  and  Eqs.  (21)  and 
(23)  that  1,  so  that  p  ~  1.  Also,  one  can  in  Eqs.  (21)  and  (24) 

generally  disregard  against  In  the  following  these  simplifications 

have  been  used.  The  integral  of  Eq.  (24)  has  been  evaluated  for  P  =  1 
and  is  plotted  in  Fig.  2.  It  will  be  referred  to  as 

Electrons,  which  are  the  primary  recipients  of  the  rf  power, 
transmit  it  partly  to  the  heavy  particles  through  elastic  and  inelastic 
collisions,  and  partly  to  the  wall  through  ambipolar  diffusion  with 
positive  ions  by  releasing  their  kinetic  energy  and  the  heat  of  recombina¬ 
tion. 

We  thus  have 

P  P  P 

r  =  <-r)  +  (?)  (25: 

gas  wall 

Disregarding  the  much  less  frequent  inelastic  collisions  and  considering, 
besides  T^,  also  the  gas  temperature  T  as  constant,  we  may  write 


(^)  =  27rR2  -1  i  k  (T  -  T) 

gas 


/*  n 

V  ‘ 


where  ^  is  the  fraction  of  its  energy  an  average  electron  loses  per 
elastic  collision  and  is  approximately  =  2m^/M,  (M  =  mass  of  heavy 
particle).  Using  the  distribution  n^(p)  of  Eq.  (11)  as  well  as  the  identity 
(17),  Eq.  (25)  can  be  written 


(g.)  =  3  7rR2  f  k 


The  portion  of  power  delivered  directly  to  the  wall  by  the  electrons  is 


11 


or  when  (dnQ/dp)^^  is  eliminated  by  the  differentiated  Eq.  (11)  and 
Eq.  (12) 

(f-)  =  ZTTD^  n  (eV  +  4-^T  )  (q  +  1)XJ,(X)  (2 

^  wall  a  o  1  i  e  i 


It  is  seen  from  Eqs.  (24),  (27)  and  (29)  that  the  total  power  as  well 
as  the  gas  and  wall  portions  are  proportional  to  ^eo*  when  these 

equations  are  substituted  into  Eq.  (24),  n^^  cancels  out.  This  permits 
us  to  evaluate  the  induced  electric  field  at  the  wall,  E^^,  as  a  function  of 
and  q  for  particular  metal-gas  mixtures  and  known  values  of  R  and  T. 


Carrying  out  the  substitution  with  the  simplified  form  of  Eq.  (24) 
and  with 

^  q+ivo3 


po  dP«l2(q) 


(3 


yields 


(q  +  i)XJ,(x)l  (3 


c)  Solution  of  System  of  Equations 

Gfeice  E^  has  been  determined  from  Eq.  (31),  the  further  pro¬ 
cedure  is  as  follows: 

E-  is  fed  back  into  Eq.  (24)  to  obtain  the  ratio  I</P 

R  O 

as  a  function  of  and  q.  To  eliminate  we 
combine  Eqs.  (12),  (14)  and  (15),  remembering  that  ™  ^gi  " 
obtain 

„  .  ^  iSJtL 

«o  J,(X)  k.(T,) 


(32 


Eq.  (32)  in  combination  with  Eq.  (24)  thus  enables  us  to  plot  P/L  as  a  function 
of  q  and  Tg.  For  a  given  value  of  P/L  we  then  obtain  a  multitude  of  couples 
qj,  T.  that  all  are  possible  solutions.  The  particular  couple  that  applies  to  a 
specific  case  is  the  one  that  also  satisfies  Eq.  (18)  which,  with  the  use  of 
Eqs.  (12),  (14),  (16)  and  (32),  can  be  written  as 


XDa(Te)  + 


4qij(q)D(T)  <| 


(33) 


Once  a  unique  value  for  q  has  been  determined,  one  has  the  complete  informa¬ 
tion  for  calculating  the  distributions  of  n  and  n^  from  Eqs.  (9)  and  (11),  respectively, 
with  aid  of  Eqs.  (12)  and  (14). 

In  the  next  chapter  application  of  the  method  to  a  particular  metal- gas 
mixture  and  specific  input  conditions  will  be  illustrated. 
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in.  APPLICATION  TO  Rb-Xe  SYSTEM 

a)  Collision  Frequency  and  Diffusion  Coefficients 

The  analysis  developed  in  the  preceding  chapter  has  been 
applied  to  an  rf  inductive  glow  discharge  in  Xe  seeded  with  Rb  which, 
as  mentioned,  is  of  interest  as  a  radiation  source  for  atomic  frequency 
standards  and  other  optically  ptmiped  devices.  For  the  particular  case 
investigated,  the  Xe  fill  pressure  is  Z  Torr  which, according  to  Ref  [21], 
corresponds  to  an  electron  collision  frequency  ~  9. 1  x  10^  sec. 

A  gas  temperature  T  =  400  K  has  been  assumed.  The  diffusion  coefficient 
of  Rb  atoms  in  Xe  has  been  calculated  from  the  well-known  expression 

D  =  (34V 

where  according  to  Me  Daniel  [25]  the  mean  free  path  of  the  Rb  atoms 
is  given  by 

/=|2^ffnd2  +  ffn  [1/2  (d +d'l^  (1  +  ^  (35) 

d  means  the  atomic  diameter  and  primed  symbols  refer  to  Xe.  Since 
n  is  several  orders  of  magnitude  lower  than  n',  the  first  term  in  Eq.  (35) 
may  be  disregarded.  With  d  ~  0.504  nm  one  obtains 

i  =  1.53  X  10"^  cm  (36) 

The  mean  velocity  of  the  Rb  atoms  follows  from 

V  =  (’^^^ )  =  3. 14  X  lO^cmsec  ^  (^"^) 

Inserting  (36)  and  (37)  into  Eq.  (34)  yields 

D  =  16.0  cm^  sec" ^  (38) 
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(39) 


The  ambipolar  diffusion  coefficient,  defined  by 

u  D.  +U.D 

D  =  -  ■?  -U.  Zf-g.  . 

»  Me  +  Mi 

can  be  approximately  calculated  from 


Da  « 


MjPe 

Me 


(40) 


A  value  for  4^,  the  mobility  of  Rb  ions  in  Xe,  has  been  obtained 
from  Loeb  [26].  At  T  =  293K  and  p  3  760  Torr  the  value  is  1.  IZcm^x 
(V  sec)~^.  Correcting  for  the  effects  of  p  and  T  on  particle  density, 
we  obtain  for  the  present  case 


=  58 1  cm^  (V  sec) 
D 

Since  — S.  =  kT. 

Me  * 


(41) 

(42) 


we  obtain  from  (40)  with  (41)  and  (42) 

=  5.01  X  10"^T^ 


(43) 


b)  Ionization  Rate  Coefficient  of  Rb  Atoms 

Tlxe  dependence  of  the  ionization  rate  coefficient  k^^  upon 
the  electron  temperature  has  been  derived  from  measurements  of  the 
ionization  cross  section  Qgi  with  low-energy  electron  beams.* 


Attempts  have  also  been  made  to  derive  values  for  kgi  more  dlrectlv 
from  the  measurements  by  Solov'ev  et  al.  [30-34]  in'^b  glow  discharges. 
These  values  are  by  up  to  six  orders  of  magnitude  below  those  from  the 
above  method  and,  in  the  present  application,  would  involve  unrealisti¬ 
cally  high  values  of  T  e«  Since  the  authors  claim  to  be  able  to  separate 
,  single-step  from  multiple-step  ionization  coefficients,  their  work  will  be 
restudied  for  the  development  of  a  three-level  discharge  model. 


I  < 
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The  qualitative  measurements  by  Tate  and  Smith  [27]  which  have  been 
calibrated  by  Brink  [28]  and  are  for  Rb  in  good  agreement  with  the  more 
recent  data  by  Me  Farland  and  Kinney  [29],  were  used. 

k^x  is  obtained  by  averaging  Q  .  over  the  range  of  electron  velocities 


where  is  the  velocity  of  an  individual  electron  and  €  is  its  kinetic 
energy.  The  averaging  is  carried  out  by  integration  over  the  electron 
energies  from  ionization  £  =  eV^  to  infinity 


where  j ^  is  the  electron  energy  distribution  function. 
is  approximated  by  the  slope  of  the  experimental  curve  at  the 
origin  «  =  <  ^ 


( 6-ej) 


as  it  is  common  practice  (see  f.  i.  von  Engel  [35]).  v^  is  given  by 


The  electron  energy  distribution  is  assumed  to  be  Maxwellian  which 
according  to  earlier  calculations  by  Margenau  [36]  represents  a  closer 
approximation  to  actual  conditions  in  rf  discharges  than  in  dc 
discharges. 


In  this  case 


where  fin=  k  T 

***  e 


Introducing  Eqs.  (46)  to  (48)  into  Eq.  (45)  then  yields 


€  *1 
udiere  X  ■  ■  —  and  X.  *  — r —  • 

« m  1  '  m 

The  integral  in  Eq.  (50)  yields 

-X. 

(X.  +  2)  e  ‘ 

From  the  experimental  data  [27  -  29]  we  have  obtained 

.  3.76x10-“^ 


(48) 


(52) 


with  the  numerical  values,  (51),  and 


X.  = 


e  Vi 


I  k  T, 


,  we  obtain  from  Eq.  (50) 


e  V 

k-.  =  2.  015  X  -z^e'iTr” 

*  I  kTe  / 


(53) 


Eq.  (53)  has  been  evaluated  with  (Rb)  =  4.  18V  and  is  plotted  in 

Fig.  3.  Also  shown  is  the  ratio  k  ./O  . 

gx  a 
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IV.  RESULTS 


a)  Data  for  Profile  Parameters 

For  calculating  the  radial  distributions  of  n  and  n^  from 
Eqs.  (9)  and  (11)  for  a  given  radius  R,  we  have  to  determine,  in 
addition  to  the  common  profile  parameter  q,  the  values  of  n^^  and  n^^. 
This  is  facilitated  by  the  charted  material  of  this  section  for  a  par¬ 
ticular  radius  R  =  0.  4  cm.  The  quantities  Ep^  and  Te.^hich  enter 
into  the.  relationships  between  dependent  and  independent  variables 
and  characterize  the  discharge  conditions,  will  also  be  presented 
graphically. 

In  Fig.  4  is  shown  the  relationship  between  Ej^,  T^.and  q  as 
obtained  from  Eq.  (30).  It  is  seen  that  T^  increases  with  Ej^,  but 
decreases  with  increasing  values  of  q.  From  Eq.  (23)  it  is  evident 
that  measurement  of  E^  can  serve  as  a  substitution  for  the  measure¬ 
ment  of  P/L. 

In  Fig.  5  P/L  is  plotted  vs.  q  with  Tg  as  a  parameter.  For 
constant  Tg,  q  increases  with  P/L,  which  indicates  that  the  profiles 
of  both  n  and  n^  will  tend  towards  a  more  rectangular  shape  as  the 
power  input  is  increased.  This  tendency  is  the  stronger,  the  higher  T^. 

The  P/L  (q,  T^  (relationship  of  Fig.  5  has  been  combined  with 
the  n^  (q;Tg)  relationship  of  Eq.  (32)  to  eliminate  Tg.  The  resultant 
n^(q;P/L)  relationship  is  plotted  in  Fig.  6.  Each  point  in  this  graph 
represents  a  solution  of  the  electron  balance  equation  for  some  particular 
input  conditions  P/L;n^,  and  it  permits  one  to  read  off  immediately  the 
appropriate  value  of  q.  The  curves  show  that  if  P/L  is  low  or  if 
n^  — n(Xe)  (s»6.6  x  10  cm”  )>q  — »0,  that  is,  distortion  of  the  profiles 
is  small.  QgCi^)  then  follows  the  regular  Bessel  function  and  n(r)  is 
uniform.  Vice  versa,  as  P/L  is  increased  or  n^  is  reduced,  the  value 
of  q  increases. 
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4.  Electron  Temperature  vs  Induced  Electric  Field  at  Tube  Wall 
for  rf  Induction  Discharge  In  Rb-Xe  Mixture 


Figure  6.  Dependence  of  Profile  Paraneter  on  Mean  Total  Rb-Partlcle 
Concentration  for  Various  Power  Inputs  per  Unit  Column 
Length  into  rf  Induction  Discharge  in  Rb-Xe  Mixture 
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The  values  of  n_  and  have  likewise  been  evaluated  as 

K  ®Q 

functions  of  P/L  and  n^  and  are  plotted  in  Figs.  7  and  8,  respectively. 

It  is  seen  from  Fig.  7  that  increases  with  both  n^.  and  P/L.  At 
n^  values  less  than  10^^/cm^  the  increase  is  approximately  pro¬ 
portional  to  as  n^  becomes  larger  the  exponent  increases, 

and  for  n^  larger  than  10^^/cm^  it  is  about  1.  The  effect  of  P/L 
on  n^^  decreases  with  increasing  n^ .  This  behavior  can  be  readily 
understood  from  the  consideration  that  the  transport  of  Rb-particles 
to  the  wall  and  the  buildup  of  concentration  require  the  expenditure  of 
energy.  The  higher  the  energy  available  per  particle,  the  higher  the 
buildup  can  go  over  the  average  concentration.  When  comparing  n^^ 
with  n^,  it  must  be  remembered  that  the  former  refers  to  atoms  only, 
whereas  the  latter  also  comprises  ions.  Thus  even  for  n^^  =  n^,  the 
concentration  of  atoms  at  the  wall  exceeds  their  average  concentration 

in  the  plasma  volume.  From  Fig.  8  it  is  seen  that  also  n^  increases 

o 

with  both  3^  and  P/L.  The  trends  here  are,  however,  opposite  to  those 
for  nj^.  The  effect  of  n^  is  largest  at  low  concentrations  and  that  of 
P/L  is  largest  at  high  concentrations  of  n^. 

In  Fig.  9  the  degree  of  ionization  as  n^/n^  is  shown  as  a  fiinction 
of  n^  and  P/L.  Within  the  range  of  parameters  considered  it  can  vary 
between  about  1%  and  35%.  As  one  would  expect,  a  increases  with  P/L 
but  decreases  with  n^.  The  intersection  of  the  curves  for  P/L  =  3  W/cm 
and  10  W/cm  is  believed  to  be  real.  To  the  higher  P/L  value  corre¬ 
sponds  a  higher  value  of  q  (Fig.  6)  and  thus  a  spreading  of  the  n^' 
distribution  closer  to  the  wall  where  the  induced  electric  field  is  higher. 
This  causes  a  higher  local  power  dissipation  (~  E^),  so  that  for  the  same 
average  ionization  one  can  have  different  power  inputs.  If  unrestricted, 
the  discharge  would  operate  at  the  highest  input,  but  in  practice  this 
negative  characteristic  of  the  discharge  is  usually  compensated  by  the 
positive  characteristic  of  the  rf  power  supply  which  forces  stable 
operation  at  some  intermediate  power  level. 
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Figure  8 


Peek  Electron  Coneentretion  at  Axis  vs  Hean  Total  Rb-Particle 
Concentration  for  Rb-Xe  rf  Discharge  Coluan 
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Figure  9 


Degree  of  Ionization  of  Rb  Particles  vs  Total  Mean  Rb-Partlcle 
Concentration  for  Inductive  rf  Discharge  In  Rb-Xe  Mixture 


28 


I 


b)  Examples  for  the  Radial  Distributions  of  Electrons  and  Rb  Atoms 

With  the  information  presented  in  Figs.  6-8  we  can  now 

calculate  absolute  radial  distributions  of  n  and  n^  from  Eqs.  (9)  and 

(11)  (with  (12))  for  arbitrary  values  of  n^  and  P/L  and  a  fixed  radius 

R  =  0.4  cm.  Such  calculations  have  been  carried  out  for  ^o  examples. 

In  the  first  one  n^  =  10*’cm  and  P/L  =  0.3  W/cm  were  chosen. 

This  yields  from  the  above  figures  q  =  2.  5,  =  2.  35  x  10^'*cm-3  and 

n^  =  1.70  X  10^ ^cm”^.  We  further  obtain  from  Figs.  5  and  4,  T,--  7500  K 
o  c 

and  1.75  V/cm.  The  corresponding  distributions  are  plotted  in 

Fig.  10.  Note  the  difference  in  scaJLes  between  n^  and  n.  It  is  seen  that 
n^  is  xiearly  constant  over  the  first  2  mm  of  tlie  radius  and  then  drops  with 
increasing  slope  towards  the  zero  value  at  the  wall.  The  distribution  of 
n  shows  the  inverse  behaviour  rising  from  a  small  fraction,  to  the  wall  value 
in  the  outer  half  of  the  radius . 

In  the  second  example  P/L  has  been  increased  by  a  factor  ten 
to  3  W/cm  for  the  same  n^  =  10^^/cm^.  This  yields  q  s>  13,  n^  =  9 
X  10^^/cm^,  n^  s4xl0^^/cm^,  T-~10^K  and  E^-^  3  V/cm.  The 
profiles  obtained  for  these  more  extreme  conditions  are  shown  in  Fig.  11. 
Compared  to  the  first  example  the  regions  of  uniform  distributions  of  n^  and  n 
have  been  extended  and  now  encompass  75%  of  the  tube  radius,  so  that 
practically  all  Rb  atoms  are  concentrated  within  one  mm  from  th^  wall. 
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Figure  10.  Radial  Dlatrlbuttona  of  Rb-Atoa  Denalty  n  and  Electron 

Density  n  In  rf  Induction-Excited  Rb~Xe  Discharge  Colunn 
P/L  -  0.3  W/c«;  (Rb)  -  10^^/cm^ 
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Figure  11. 


Badlel  Dlstrlbutloae  of  Rb^'Atom  Density  n  and  Electron 
Density  n.  in  rf  Induetlon-bcited  Rb^Xe  Discharge  Coluan 
P/L  -  3  W/ci«  \  (Rb)  -  lO^^/cm’ 
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DISCUSSION 


There  are  at  present  no  experimentally  determined  particle  profiles 
or  related  data  from  induction  discharges  with  which  the  material  of  Figs.  6-11 
can  be  compared.  The  closest  data  for  comparison  are  the  measurements  by 
Dutch  workers  in  the  positive  column  of  Cs  -  Ar  dc  discharges.  The  shape 
of  the  n  -  profiles  obtained  by  Waszink  and  Polman  [4]  is  similar  to  those  in 
Figs.  10  and  11.  Also,  the  Cs-atom  density  profiles  by  Bleekrode  and 
V. d.  Laarse  and  by  van  Tongeren  [S,  6^  agree,  for  the  higher  discharge 
currents,  qualitatively  with  the  Rb  profiles  in  the  above  figures.  The  trend 
towards  a  more  rectangular  profile  shape  with  increasing  P/L,  evident  from 
Figs.  10  and  11,  is  observed  in  the  above  experiments  as  a  function  of  discharge 
current. 

With  respect  to  Rb  discharges  in  Ne,  Ar,  or  Kr,  Tako  et  al.  [18]  and  Hirano 
(193  concluded  from  their  observations  on  spherical  lamps  that  emitting  as  well 
as  absorbing  Rb  atoms  are  concentrated  near  the  walls. 

Regarding  the  electron  temperatures,  preliminary  experiments  at 
The  Aerospace  Corporation  with  Rb-Xe  induction  discharges  of  the  same 
radius  as  used  in  the  calculations  indicate  that  the  calculated  values 
may  be  too  high  by  a  factor  of  two  or  three.  Several  features  in  the  model 
may  contribute  to  such  a  discrepancy.  The  major  cause  is  probably  an 
underestimation  of  the  electron  production  rate  by  disregarding  multi- 
step  ionization,  especially  two-etep  ionization  via  the  first  excited  state  [7,  S.J 
Such  multiple-step  processes  become  important  at  higher  ionization  levels. 

They  require  adding  terms  proportional  to  n  and  higher  powers  in  n  to 
the  electron  balance  equation  (52)  and  thus  make  this  equation  nonlinear. 

In  addition  to  the  underestimation  of  production  rates  it  is  also 
likely  that  the  model  in  its  present  form  overestimates  electron  -  Rb  ion 
loss  rates.  The  assumed  distribution  of  Rb  atoms  according  to  Eq.  (9) 
involves  n^  -  0  for  all  values  of  q>0,  so  that  there  is  no  electron  production 
at  the  axis. 
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This  condition  will  be  approached  for  high  values  of  P/L,  but  at  low  inputs 
we  must  expect  n^  to  be  nonnegligible.  Electrons  produced  at  the  axis  will 
have  the  longest  distance  to  travel  to  the  wall;  they  will  thus  in  the  average 
have  the  longest  life  in  the  discharge  and  need  to  be  replaced  at  a  slower  rate. 
Excltiding  electron  production  from  the  central  region  and  shifting  it  closer  to 
the  wall  will  require  higher  production  rates  which  are  associated  with  higher 
values  of  and/or  n^. 

Allowing  for  a  finite' concentration  of  Rb  atoms  at  the  axis  accord¬ 
ing  to 

*  —  X 

(9a) 

i 

would  make  the  assumptions  in  the  model  more  consistent  for  low  P/L  and 
large  n^. '  Since  the  addition  of  n^^  0  would  leave  the  production  rates  linear 
in  n  under  otherwise  unchanged  conditions  (no  consideration  of  multistep 
ionization), the  solution  for  the  space-dependent  term  in  Eq.  (9a),  which  is 
Eq.  (11),  could  be  superimposed  on  the  solution  for  the  constant  term  which 
is  the  ordinary  Bessel  function.  The  resulting  increase  in  complexity  of  the 
analysis  would  be  moderate. 

The  electron  temperature  does  not  appear  explicitly  in  the  data  of 
Figs.  6-11.  It  will  affect  them  mainly  through  its  influence  on  (£^q.  43). 
Having  this  value  too  high  will,according  to  Eq.  (15),  yield  too-high  levels 
of  2  relative  to  n  . 

Aside  from  these  apparent  deficiencies  of  the  present  model,  it  is 
believed  to  give  a  qualitatively  correct  picture  of  the  conditions  in  an 
alkali  metal  -  noble  gas  discharge  and  of  the  relationships  involved. 

A  conclusion  that  can  be  drawn  from  the  results  with  respect  to  the 
retention  of  Rb  particles  at  the  wall,  mentioned  in  the  introduction,  is 
that  lowering  the  total  Rb  particle  concentration  n^  (by  lowering  the  vapor 
pressure)  does  not  cause  a  proportionate  decrease  in  the  wall  concentration 
np,  since  the  remaining  particles  are  driven  that  much  harder  to  the  'wall. 
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LA80RAT0RY  OPERATIONS 


The  L«boratory  Operetions  of  The  Aerotpece  Corporation  la  conducting  exper- 
ioental  and  theoretical  Inveatlgatlona  neceaaary  for  the  evaluation  and  applica* 
cion  of  scientific  advances  to  now  ttllitary  apace  ayateaa.  Veraatlllty  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  peraonnel  Ir. 
dealing  with  the  laany  probleas  encountered  in  the  naclon’a  rapidly  developing 
space  systeas.  Expertise  in  the  latest  scientific  developaents  is  vital  Co  the 
accoapllshaent  of  casks  related  to  these  probleas.  The  laboratories  chat  con¬ 
tribute  to  this  research  are: 

Aerophyslcs  Laboratory;  Launch  vehicle  and  reentry  aerodynaalcs  and  heat 
transfer,  propulsion  chealstry  and  fluid  aechanlcs,  structural  aechanlcs,  flight 
dynaaics;  high-teaperature  thenaoaechanlcs ,  gas  kinetics  and  radiation;  research 
in  environmental  chemistry  and  contaalnation;  cw  and  pulsed  chemical  laser 
development  including  chealcal  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  ataospheric  propagation,  laser  effects  and  countermeasures* 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattetlng,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bioenvironmental  research  and 
monitoring. 

Electronics  Research  Laboratory;  Microelectronics,  CaAs  low-noise  and 
power  devices,  eemiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  lidar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  Imaging;  millimeter-wave  and  microwave  technology* 

Information  Sciences  Research  Office:  Program  verification,  program  trana- 
latlo*?ri  performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  Intelligence, 
and  microelectronics  applications* 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
space  environment;  materials  performance  in  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  in  enemy-induced  environments* 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  phyaics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  airglow;  magnetoapheric  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  ^n  space  systems* 


